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Abstract 

PURPOSE: To provide a method for screening continuous color tone images. 
CONSTITUTION: A non-processed image pixel is selected according to a random spatial filling 
two-dimensional curve, a regenerative value H(ij) to be used for recording the image pixel on a 
recording medium 80 is determined from a color tone value P(i j) of the non-processed image 
pixel, an error value E(iJ) is calculated based on the difference between the color tone value of the 
non- processed image pixel and the regenerative value, and the error value is added to the color 
tone value of the non-processed image pixel by a block 70. Next, the color tone value is replaced 
with the provided total value or the color tone value of each non-processed image pixel to 
distribute the error value is replaced with the total of the color tone value of the non-recessed 
image pixel and the part of error so that the error value can be distributed over more than two of 
the non-processed image pixels. These steps are repeatedly processed until all the image pixels 
are processed by a counter 8. 
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(57)Abstract: 

PURPOSE: To provide a method for screening continuous color 
tone images. 

CONSTITUTION: A non-processed image pixel is selected 
according to a random spatial filling two-dimensional curve, a 
regenerative value H(i,j) to be used for recording the image 
pixel on a recording medium 80 is determined from a color tone 
value P(i j) of the non-processed image pixel, an error value E 
(i j) is calculated based on the difference between the color 
tone value of the non- processed image pixel and the 
regenerative value, and the error value is added to the color 
tone value of the non-processed image pixel by a block 70. 
Next, the color tone value is replaced with the provided total 
value or the color tone value of each non-processed image 
pixel to distribute the error value is replaced with the total of the 
color tone value of the non-recessed image pixel and the part of error so that the error value can be 
distributed over more than two of the non-processed image pixels. These steps are repeatedly 
processed until all the image pixels are processed by a counter 8. 
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typedef struct { 
int i; 
int j; 

> Index; 

typedef struct { 
rndex ptl; 
Index pt2; 
Index ptS; 
Index pt4; 

> Hi lbert__Ele»; 

typedef struct. { 

char type [80]; /* type description */ 

char dascr [80}; /* variable description*/ 

int nry /* number of rows */ 

int nc; /* number of columns */ 

int min; /* min value range */ 

int max; /* max value range */ 

^ n t **el«n; /* elements of matrix */ 

) Itile; 

store_hilbert_elem(it, p) 
Itile *it; 
Hilbert_JBlem *p; 

{ 

static int n=0; 

it->elem[n] [0J -p->ptl-i; it->elem[n] [1] »p->ptl.j; n++; 
it->elexn[nj [0] - p->pt2.i; it->elem[n] [1J - p->pt2.}; n++; 
it->elem[n] [0] - p->pt3 r i; it->elem[n] [1] = p->pt3.;J; n++; 
it->elem[nH0] - p->pt4.i; it->elem[n] [1] = p->pt4.j; n++; 
) 



C3t2] 
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/*** Initiation for Recursive Calculation of Hilbert Scan ***/ 



hilbert_Lnitiation(size,p) 
lnt siza; Hilbert_Elera *p; 

< 

p->ptl.i - l*siz*/4; p->ptl.J » l*size/4; 
p->pt2.i - l*3lze/4; p->pt2. j . - 3*size/4/ 
p->pt3.i - 3*size/4; p->pt3 . j » 3*size/4; 
p->pt4.i » 3*size/4; p->pt4.j = l*size/4; 
} 

/*** Recursive Module to Calculate Hilbert Scan ***/ 
hilbert_j>ropagation (lt,p) 

Itile *it; 

Hilbert_Elem *p; 

{ 

int il r jl,i2, j2, 14, j4, 1, li, Ij, si, s j, orientation, length; 
Hilberc_Elem pl,p2 r p3,p4; 



il 




p->ptl.i/ 


jl 




p->ptl.j/ 


12 




p— >pt2.i; 


32 




p->pt2. j; 


14 




p— >pt4 . i ; 




a 


p->pt4.j; 



li - (i4 - il)/2.0,- 
lj - <j4 - jl)/2.0; 

orientation- (i4-il> * ( j2-jl) - ( j4-jl) * (i2-il) ; 
1 « (int) sqrt ( {double) li*li + 1 j*l j] / 
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g 1 

if (orientation < 0) 

( 

pl.ptl.i - p->ptl.i; 
pl.ptl.i — <li+lj)/2.0; 
pl.pt2.i » pl.ptl.l+li; 
pl.pt3.i « pl.pt2.i+lj; 
pl.pt4.i - pl.pt3.i-li; 

p2.ptl.i = pl.pt4.i+lj; 
p2.pt2.i - p2.ptl.i+lj; 
p2.pt3.i - p2.pt2«i+li; 
p2.pt4.i - p2.pt3.i-lj/ 

P 3.ptl.i - p2.pt4.i+li; 
p3.pt2.i « p3.ptl.i+lj; 
p3,pt3.i - p3.pt2.i+li; 
p3.pt4.i = p3.pt3.i-lj; 

p4.ptl.i =» p3.pt4.i-lj; 
p4,pt2.i - p4.ptl.i-li; 
p4.pt3.i - p4.pt2.i-lj; 
p4.pt 4. i « p4.pt3.i+li; 
) 

else 
< 

pl.ptl.i - p->ptl.i; 
pl.ptl.i — ' <li+lj)/2.0; 
pl.pt2.i = pl.ptl.i+li; 
pl.pt3.i - pl.pt2.i-lj/ 
pl.pt4.i = pl.pt3.i-li/ 

p2.ptl.i - pl.pt4.i-lj; 
p2.pt2.i « p2.ptl.i-lj; 
p2.pt 3.1 « p2.pt2.i+li; 
p2.pt4.i = p2.pt3.i+lj; 
[14] 40 




ftmy-G -7 0 14 4 

20 

*<o 3 ) 



pl.ptl.j « p->ptl.j; 
pl.ptl.j — (li+lj)/2.0; 
pl.pt2.j - pl.ptl.j+lj; 
pl.pt3.j - pl.pt2.j-li; 
pl.pt4. j - pl.pt3. j-lj; 

p2.ptl.j « pl.pt4.j-li; 

p2.pt2.j - p2.ptl.j-li; 

p2.pt3.j = p2.pt2.j+lj/ 

p2.pt4.j - p2.pt3.j+li; 

p3.ptl.j - p2.pt 4. j+lj; 
p3.pt2.j « p3.ptl.j-li; 
p3.pt3.j « p3.pt2.j+lj; 
p3.pt4,j « p3.pt3. j+li; 

p4.ptl.j - p3.pt4.j+li/ 
p4.pt2.j - p4.ptl.j-lj; 
p4.pt3.j = p4.pt2.j+li; 
p4.pt4.j » p4.pt3. j+lj/ 



pl.ptl.j « p->ptl. j; 
pl.ptl.j — (li+lj)/2.0; 
pl.pt2.j = pl.ptl. j+lj; 
pl.pt'3-j = pl.pt2. j+li; 
pl.pt4.j - P L.pt3. j-lj; 

p2.ptl.j « pl.pt 4. j+li; 
p2.pt2.j « p2.ptl.j+li; 
P 2.pt3.j - P 2.pt2. j+lj; 
p2.pt4.j = p2.pt3. j-li; 
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m 1 



(*<D If) 



p3.ptl.i = 
p3.pt2.i - 
p3.pt3.i - 
p3.pt4.i " 



p2.pt4.i+li; 
p3.ptl.i-lj; 
p3.pfc2.i+li; 
p3.pt3.i4lj; 



p3.ptl.j - p2.pt4. j+lj; 

p3.pt2.j » p3.ptl.j+li; 

p3.pt3.j = p3.pt2.j+lj; 

p3.pt4, j - p3.pt3.j-li; 



p4.ptl.i - p3.pt4.i+lj; 
p4.pt2.i — p4.ptl.i— li; 
p4.pt3.i « p4.pt2.i+lj; 
p4.pt4.i *■ p4.pt3.i+li; 



p4.ptl.j •* p3,pt4 . j-li; 
p4.pt2.j = p4.ptl.j— Ij; 
p4.pt3.j « p4.pt2. j-li; 
p4.pt4.j =» p4.pt3.j+lj; 



if (1 > 1.0) 

{ hilbert_propagafcion (it, &pl) ; 
hilbert_jpropagation (it , £p2) ; 
hilbert_jpropagation (it, 6p3) ; 
hilbert_propagat ion ( i t , &p4 ) / } 
else /*" termination */ 
{ stojre_hilbert_elem (it , 6pl J ; 
store^hilbert^ elem< it , &p2) ; 
atore_hilbei-t_elem(it, &p3) ; 
store_Jiilbert_elero(it, fcpl) ; 
retum; } 



main () 
{ 

char name_j)ath [ 32 ] ; 
int. size; 
FILE *fp; 
Itile it; 



[&5] 
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Hilbert_Elera p; 

printf ("enter name path under which the Hilbert path will be stored: 



acanf (°%s" / name_path) ; 

tp « £open(name_path,"w") ; 

printf Center size of square path (in pixels, must be power of 2! !): 
"># 

scanf ( n %d n / Ssi«e) ; 
size =° 32; 

alloc_itile (size*size, 2, Sit) ; 
atrcpy (it .descr, w hilbert_curve") ; 
it.ru: - slze*3ize; 
it.nc « 2; 
it .ruin = 0; 
it. max — size? 

hilbert_initiation(size, tp) ; 
hilbert_jpropagation(&it, &p) ; 
write__itile (fp,£it) ; 
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permit^ 2D(seed f n, a) 
int aeed r »,**a; 

{ 

int i,*b,ct2],d[21; 

b « (int *) ivector <n*n) ; 

ran_perturb <sead f n*n,b) ; 

/* replaces the n x n elements in vector b by a random permutation*/ 
el01-*c[l]*"n/ 
for (i«0; i<n*n; i+f) 
i 

calc_index__ £ rom_lin_ addr (2, c,b [i] , d) ; 
/* transform the linear address b[i] into a coordinate pair in 
vector d */ 

a[d[0]] IdCin - i; 

> 

free_i vector (b) ; 
) 

/*** RECURSIVE CAXCU1ATI0N OF 2D ORDER ***/ 
recura^order^calc <sd, lv, tp, nb, ib, jb, od) 

int *sd,lv,*tp,nb,ib, jb; 

Itile *od; 

{ 

int i, j,**ma,sz,ba; 
az. - tp[lv] ; 

ma - (int **) imatrix(tp[lv] ,tp(lv] ) ; 
permut_ 2D Cad, sr,ma) ; 
for ( i-0 , ba-1 ; i<lv; i++) 

ba *- tpCi] ; 
if(lv — 05 
( 

fox(i~0;i<sz;i++) 
for ( j-0; j<sx; 



od->Gl«nCnb4-ma[i] [j]*ba*ba] £0] - ib+i*ba; 

od->e^emtnb+inaCi] I j]»ba*ba] [1] = jb+j*ba; 
) 

return; 
) 
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g 2 -2J 

£rea — imatrix (»fi f mo) ; 
> 

main () 
I 

char namejpath[32] ; 

int n / 3eod,l«vel / topol[S] f **order; 

int site; 
FILE *fp; 
Itile it; 

print* ("enter name of path: 

scan£( ,T *3",namejpatli) ; 

£p - f open (namejpath, w w M ) ; 

/* -.i..- la the of matrix over which error propagation will 

take place */ 

/t'thls^trix will be "level- time* recurai.ely aubdivided into 

sub square 3 */ 

level - 4 ; . 
/* -topol- deacribes the aubaequent aize of tfceae autanitricea / 
topol[0]-2; topoltll-2, fpol[2]-2» *opol(3]-2» topol[4]=2,- 

Blloc_itila(sis«*aize,2,6lt) ; 

atrepy (it .desor, "cr^ath") t 

it.nr — sire*»ize; 

it .nc - 2; 

it .min — 0; 

it .max — aize; 

seed - -1; 

recur 3 _ordet_c»lc (5 seed, level . topol, 0 , 0 , 0 , Sit) i 

write_itile (fp } Sit) ; 

) 

[021 S2 F SSiSfc"* HWWtWanTM. 
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1 3 






2 0 






3 0 






3 1 






4 0 


a — y h—Xkg© 




4 1 




5 vn v >; 


3 3, 


4 3 ki/7M/vX^ 


8 mMti t y>& 


6 3 


li^hi^X* 


9 tyaD4' umMmm 


6 0 




1 0 LUT i0 


5 0, 


6 5, 7 0 ifcfitgB 


1 2 *tiHHBtfe£B 


8 0 





[01] [02] 
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